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ABSTRACT 

Spectral reflectance of snow under diffuse illumination is studied 
using the two-stream approximation of the radiative transfer equation. The 
scattering and absorption within the snowcover due to the randomly distrib- 
uted ice grains are characterized by the single scattering albedo and aniso- 
tropic phase function. Geometric optics calculations are used to relate the 
scattering and absorption parameters to the grain size and density of snow. 
Analytical expressions for the intensity within the snowpack and the asymp- 
totic flux extinction coefficient are also obtained. Good agreement is shown 
between the theory and available experimental data on visible and near- 
infrared reflectance, and the asymptotic flux extinction coefficient. The 
theory also may be used to explain the observed effect of aging on the 
snow reflectance. 
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TWO-STREAM THEORY OF SPECTRAL 
REFLECTANCE OF SNOW 

INTRODUCTION 

Interest in developing techniques for remote measurement of snow parameters (e.g M density, 
grain size, water content) have increased significantly in recent years. One of the more 
promising techniques to monitor these parameters is that of microwave radiometry. A recent 
study by Chang and Choudhury ( 1 978) has used a radiative transfer model to explain the 
observed brightness temperature of polar flm from the measured physical temperature, 
density and grain size. Analysis of spectral reflectance is an alternative, and complementary 
to microwave radiometry. Apart from remote sensing purposes, the study of spectral reflec- 
tance is relevent in climatology because it determines the amount of solar radiation absorbed 
by the earth’s surface. 

There have been several theoretical attempts to relate the spectral reflectance of a model 
snowcover with its physical properties (e.g., Dunkle and Bevans, 1 956; Giddings and La Cha- 
pelle, 1961 ;Bergen, 1970, 1971, 1975; Bohren and Barkstrom, 1974). These theories either 
contain parameters which cannot be obtained directly from the measured physical proper- 
ties or are valid for a restricted range within the visible spectrum. 

For visible and near-infrared radiation, the scattering by individual ice grains is such that a 
major portion of the radiation gets scattered in the forward direction (Bohren and Barkstrom, 
1974). The high reflectivity of snow is due to the multiple scattering effect within snow. 
Calculation of spectral reflectance necessarily involves a study of radiative transfer in snow. 
For this purpose, the radiative transfer equation is useful if snow is represented by ran- 
domly distributed, well separated ice grains. The snow parameters enter into the radiative 
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transfer equation through the single scattering albedo, the extinction coefficient and the 
phase function. 

The physical properties of a natural snowcover are generally not known accurately and are. at 
best, approximately homogeneous. Although an accurate solution of the radiative transfer 
equation will be useful for standard reference and comparison purposes, to date no such cal- 
culation has been performed. Based on an approximate solution of the radiative transfer 
equation, a simplistic theory of spectral reflectance is given in this paper. All parameters in 
the theory can be calculated directly from the measured physical properties of snow, and the 
theory is applicable to visible and near-infrared radiation. The exact points of departure be- 
tween the analysis presented below and the previous reported results vary. Instead of dis- 
cussing these differences, a comparison will be made to show the relationship among the re- 
sults. In addition, comparisons will be made between the theory and the observations. 

APPROXIMATE EXPRESSION FOR REFLECTANCE 

To obtain a closed form analytic expression for reflectance, the following assumptions will 
be made: 

1 . The snowcover is homogeneous. 

2. Radiative transfer within the snowcover can be studied by dividing the radiation 
field into two streams of intensities. 

3. Geometric optics calculations can be used to obtain the scattering and absorption 
parameters of individual ice grains. 

With the first assumption the presence of internal inhomogeneities is neglected. This 
assumption is particularly crucial for a natural snowpack because of metamorphism. This is 
also a common assumption in all previous models. 




The two-stream representation of the radiation field is a widely used method, and shown to 
be a reasonably accurate approximation of the radiative transfer equation (e.g., Oomoto and 
Wang, 1974;Coakley and Chylek, 1 975 ; Wiscombe and Grams, 1976; Venkatram and 
Viskanta, 1977). Although multiple scattering problems can be solved analytically with this 
approximation, the use of this approximation is a weakness of this study. A special form 
of this approximation was used previously in the calculation of snow reflectance (Dunkle 
andBevans, 1956; Bergen, 1975). 

In general form, the two-stream approximation corresponds to replacing the radiative trans- 
fer equation by the following coupled differential equations for mean intensities (I + (r) and 
I.(r)) in the forward and the backward hemispheres: 

dl + dl 

6— = “1 "0) 1 + +w0I ! = +«(1-WI + M0I. 

dr dr “ 

where w is the single scattering albedo, and r is the optical depth within the snow, given by 


r = 7 e z 

where 7 e is the extinction coefficient, and z is the depth measured from the snow surface 
within the snowcover. The coefficients 5 and 0 are parameters, the choice of which distin- 
guishes alternate forms of the two-stream approximation. Some of the choices of 6 and 0 
are as follows: 

(i) Generalized Schuster-Schwartzschild (Wiscombe and Grams, 1976) 

01 *T f d/1 /*d/i' p(/i,-/i') 

*-Jo Jo 
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(ii) Modified Schuster-Schwartzschild (Sagan and Pollack, 1967; Lyzenga, 1973) 


«2 


l 



l i t 

0 2 = - 1 1 - —J p(cos 0) cos 0 sin 6 60 

where p(*i, n') and p(cos 6) have usual meaning of the phase function (Chandrasekhar, 
1960). 


It is expected that different set of 5 and 0 will give different numerical results. The merit of 
a choice for a particular application should be decided by comparison with the exact solution 
of the radiative transfer equation. Unfortunately, as stated in the introduction, such a com- 
parison cannot be made for the present application. In the following, the method of deter- 
mining the parameters and the solution of the equations will be discussed. 

The phase function needed for the calculation of 0 can be obtained from experimental 
observations or from a scattering theory calculation. Since there are no available experi- 
mental observations, one needs to perform a theoretical calculation by specifying the shape, 
size and the refractive index of ice grains. By discussing the difficulty and the limitations 
associated with the choice of shape, a recent study on snow (Bohren and Barkstrom, 1974) 
considered spherical non-absorbing particles and used geometric optics to calculate the 
phase function. Geometric optics calculations for transparent spheres (Bohen and 
Barkstrom, 1974; Hansen and Travis, 1974) give the value of the parameter 0 2 as 0.065. 

Using the approximate relationship between the parameters 0, and 0 2 (Wiseombe and 
Grams, 1976) the value of the parameter 0, is obtained as 0.12 ± .05. The uncertainty in 0 2 
is due to the stated accuracy of the approximate relationship. Because of this uncertainty, 
only the modified two-stream approximation will be compared with observations. The single 
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scattering albedo, w, can be obtained from the same calculational method which gives the 
phase function. A simple formula which reproduces the single scattering albedo with con- 
siderable accuracy as compared with the Mie theory results is (Irvine and Pollack, 1 968): 

w= T + i expH ' 67 k * r) 

where k x is the absorption coefficient of ice (Hobbs, 1974; Irvine and Pollack, 1968), and r 
is the radius of the sphere. For non-spherica] particles, one can use this formula by modify- 
ing the numerical coefficient in the exponential (Sagan and Pollack, 1967). 

The extinction coefficient of snow y e is the product of the extinction cross-section and the 
number density of ice grains. By taking the extinction cross-section as 2rrr 2 , one obtains 



where p and p i are respectively the density of snow and ice. 

All parameters in the radiative transfer equation have now been related with the physical 
properties of the snowcover. 

If the surface of the snowpack is not ice-glazed, then for diffuse incident radiation, I 0 , the 
boundary conditions for the solution of the equations are: 

M°)-i 0 

I_(r 0 ) = RI + (r 0 ) 

where R is the reflectivity of the underlying surface and r 0 is the total optical thickness of 
the snowcover, 
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f 0 a ^ h ' 


where h is depth of the snowcovcr 

Following standard procedures (Sobolev, 1960) the solution can be written as: 


where 


C| p i - ojj ^2 r i * j 

C) r i - wl ^2 f i - oj~ | 

2 5a J 2|_ 5a J 


a 2 = 


( 1 - u>X 1 - w + 2w/J) 


5 2 


2L 


c > " 


1-col 

1 4- - 

[, + «(a-R) 2arn] 

5a 

Ra- 1 J 
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\Ra- 1/ 


1 0 e 2ar o 


[ ,+i ir][ 
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Ra- 1 


a= 1 - 


2(1 - cor 


(1 - w + 2w0)* +(1 - w) v: 


From this solution, the reflectance A(r 0 ) can be calculated as: 

I (0) 

A(r„> = — 
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(Rq- l)a ■*•(«- R)e~ 2 * r o 
(Ra- 1) + a(a- R)e" 2af ° 

For a deep snowpack (r 0 -* °°), the reflectance is given by 


A(«) = a 


The asymptotic flux extinction coefficient is obtained as: 



COMPARISON WITH PREVIOUS CALCULATIONS AN.) OBSERVATIONS 
Tne reflectance of a deep snowcover and the asymptotic tli x extinction coefficient can be 
calculated from the grain size and the density of the snowcover. Besides the assumptions 
made in obtaining the analytical results, the model of split : cal particles was used to cal- 
culate the optical properties of individual ice grains. Although the back-scattered fraction, 

0 . is shown not very sensitive to the particle geometry (Wiscombe and Grams, 1976), other 
optical parameters (i.e., the single scattering albedo and the extinction coefficient) do 
depend upon the geometry. The ice grains of a snowcover are not spherical. The shape is 
generally oblong and of variable dimensions. Although average optical properties of the 
grains, due to their random orientation with respect to the mean intensity, may correspond 
to a sphere, reservation should be exercised regarding the applicability of the theory to a real 
snowpack. Thus, the comparison between the theory and the experimental observation 
presented below should not be regarded as conclusive. 
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Bergen (1970, 1975) related semi-empirically the snowpack parameters to the reflectance, 
and to the asymptotic flux extinction coefficient derived by Dunkle and Bevans (1956). 
Using the approximate solution of the radiative transfer equation derived by var. de Huist 
(1970, a, b), and taking the model of spherical ice grains, Bohren and Barkstro/n (1974) 
related the snw parameters to the reflectance, and to the asymptotic flux extinction co- 
efficient for visible radiation. Since Bohren and Barkstrom have compared their theory 
with Bergen’s results, the following comparison will be restricted to Bohren and Barkstrom’s 
theory. The status of the experimental observations of the reflectance and the flux extinc- 
tion coefficient is disappointing either because of quantitative disagreement (Lilljequist, 
1956; Thomas, 1963) or due to incomplete specification of snow parameters (O’Brien and 
Munis, 1975). Since Bohren and Barkstrom have compared with Lilljequist’s observations, 
this comparison will be discussed first with the present theory. Figure l and Table 1 show 
this comparison for the asymptotic flux extinction, and the reflectance of a deep snow- 
cover. Although present theory appears to be little inferior to Bohren and Barkstrom's 
theory, it is difficult to vindicate either theory because the observations have not yet been 
duplicated. The discrepancy between the models is particularly noticeable in the asymptotic 
flux extinction coefficient. Whereas Bohren and Barkstrom predict lower values, the present 
theory predicts higher ones than the observations. 

Figure 2 we show the calculated snow reflectance compared with the experimental values 
for a nearly fresh snow (O’Brien and Munis, 1975). The results of the calculation are in 
good agreement with the observations. All prominent spectral structures appearing in the 
observation are well duplicated in the calculations. The quantitative agreement is 
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also good but should be treated with caution because the experimental results are 
relative to a standard (white Barium Sulfate) reflector. It is however reassuring that 
the radius values used in the calculation are in the range expected for a fresh fallen snow. 

The ice crystals of fresh fallen snow usually are of complex form and contain sharp 
corners. The shape of these crystals changes with time depending upon the vapor con* 
tent and the prevalent temperature. The process of equi*temperature metamorphism 
leads to the production of fairly uniform and well rounded grains. At the initial stage 
of this metamorphism the mean radius of the ice crystals is about 0.2 mm. The radius 
then continues to increase as the process of metamorphism advances. At a fairly 
advanced stage of this metamorphism, the radius increases to about 1.0 mm. The effect 
of other metamorphism such as the temperature-gradient metamorphism or the melt- 
freeze metamorphism, is generally to produce non-spherical and non-uniform ice crystals 
also the crystals are larger than equi-temperature metamorphism. Thus, although it is not 
unique, one should be able to use the radius of ice crystal to characterize the stage of 
metamorphism. Figure 3 illustrates the changes in the calculated snow reflectance due to 
the difference in the crystal radius. The chosen radius an* the typical snow crystal sizes 
for different stages of the equi-teinperature metamorphism. The overall shape of the 
curve does not seem to depend upon the radius of the ice crystals but the relative magni- 
tudes of various parts of the curve vary with the radius. The spectral reflectance in the 
red and near-infrared regions shows maximum sensitivity to the crystal sizes. 

CONCLUSION 

Based on the two-stream approximation of the radiative transfer equations, analytic results 
for the reflectance and the asymptotic flux extinction coefficient of a homogeneous snow- 
pack are derived. Using geometric optics calculations for spherical particles, the reflectance 
and the extinction coefficient are related to directly measurable snow parameters (density, 
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and the grain size). Good agreement was shown with the observed spectral dependence of 
the reflectance and the extinction coefficient using the measured (expected) snow param- 
eters. Further observations are highly desirable to test the accuracy and hence to refine the 
calculation. 
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Table 1 


Comparison of Observed and Calculated Values of the Reflectance 
Snow Parameters - p- 0.43 gm/cm 3 , r * 0. 1 5 mm) 


Wavelength 

(M) 

Modified 

two-stream 

Observation* 

Bohren 

and 

Barkstrom 

0.420 

0.987 

0.960 

0.978 

0.520 

0.980 

0.972 

0.967 

0.595 

0.974 

0.960 

0.954 

0.650 

0.964 

0.921 

0.940 


♦Data, UlJjequist (1956). 
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★ DATA: ULLJEQUIST (1966) 



Figure t. Comparison of Observed and Calculated Values of the Asymptotic 
I'lux 1'xtinction Coefficient. 
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